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We consider a system consisting of a particle in the harmonic approximation, having frequency 
lj, coupled to a scalar field inside a spherical reflecting cavity of diameter L. By introducing 
dressed coordinates we define dressed states which allow a non-perturbative unified description of 
the radiation process, in both cases, of a finite or an arbitrarily large cavity, for weak and strong 
coupling regimes. In the case of weak coupling, we recover from our exact expressions the well known 
decay formulas from perturbation theory. We perform a study of the energy distribution in a small 
cavity, with the initial condition that the particle is in the first excited state. In the weak coupling 
regime, we conclude for the quasi-stability of the excited particle. For instance, for a frequency 
uj of the order uj ~ 4.00 x 10 14 /s (in the visible red), starting from the initial condition that the 
particle is in the first excited level, we find that for a cavity with diameter L ~ 1.0 x 10 _6 m, the 
probability that the particle be at any time still in the first excited level, will be of the order of 97%. 
The behaviour of the system for strong coupling is rather different from its behaviour in the weak 
coupling regime. For appropriate cavity dimensions, which are of the same order of those ensuring 
stability for weak coupling, we ensure for strong coupling the complete decay of the particle to the 
QQ ' ground state in a small ellapsed time. Also we consider briefly the effects of a quartic interaction up 

^\ ' to first order in the interaction parameter A. We obtain for a large cavity an explicit A-dependent 

expression for the particle radiation process. This formula is obtained in terms of the corresponding 
exact expression for the linear case and we conclude for the enhancement of the particle decay 
induced by the quartic interaction. 
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I. INTRODUCTION 



Exact solutions of problems in theoretical physics is known to researchers since a long time ago to be a rather 
rare situation. It is a common feature to different branches of physical sciences, such as celestial mechanics, field 
theory and statistical physics, that exact solutions of coupled equations describing the physics of interacting bodies 
is a very hard problem. In statistical physics and constructive field theory, general theorems can be derived using 
cluster-like expansions and other related methods [1]. In some cases, these methods allow the rigorous construction of 
field theoretical models (see for instance [2] and other references therein), but, in spite of the rigor and in some cases 
the beauty of demonstrations, these methods do not furnish useful tools for calculations of a predictive character. 
Actually, apart from computer calculations in lattice field theory, the only available method to solve this kind of 
problems, except for a few special cases, is perturbation theory. In modern physics, a prototype situation for instance 
in abelian gauge theories, is a system composed of a charged body described by a matter field interacting with a 
neutral (gauge) field through some (in general non linear) coupling characterized by some parameter g, usually named 
the coupling constant or the charge of the body. The perturbative solution to this situation is obtained by means of 
the introduction of bare, non interacting matter and gauge fields, to which are associated bare quanta, the interaction 
being introduced order by order in powers of the coupling constant in the perturbative expansion for the observables. 
This method gives remarkably accurate results in Quantum Electrodynamics and in Weak interactions. In high energy 
physics, asymptotic freedom allows to apply Quantum Chromodynamics in its perturbative form and very important 
results have been obtained in this way in the last decades. 

In spite of its wide applicability, there are situations where the use of perturbation theory is not possible, as in 
the low energy domain of Quantum Chromodynamics, where confinement of quarks and gluons is believed to take 
place. In this particular situation, no analytical approach in the context of Quantum Field Theory is available up to 
the present moment. Also there are situations in the scope of Quantum Electrodynamics, in the domain of Atomic 
Physics, Cavity Quantum Electrodynamics and Quantum Optics, where perturbation methods are of little usefulness. 
The theoretical understanding of these effects on perturbative grounds requires the calculation of very high-order 
terms in perturbation series, what makes standard Feynman diagrams technique practically unreliable in those cases 
[3]. In this article we study situations where these mathematical difficulties can be circumvected. We consider systems 
that under certain conditions may be approximated by the system composed of a harmonic oscillator coupled linearly 
to the modes u>i of a scalar field trough some effective coupling constant g (2), the whole system being confined in 
a cavity S. We present an unified formalism for the truly confined system in a small cavity and for the system in 



free space, understood as the limit of a very large cavity. Linear approximations are used in several contexts, as in 
the general QED linear response theory, where the electric dipole interaction gives the leading contribution to the 
radiation process [4,6] In cavity QED, in particular in the theoretical investigation of higher-generation Schrodingcr 
cat-states in high-Q cavities, a linear model is employed [7]. Also, approaches of this type have been used in condensed 
matter physics for studies of Brownian motion [8,9] and in quantum optics to study decoherence, by assuming a linear 
coupling between a cavity harmonic mode and a thermal bath of oscillators at zero temperature [10,11]. In another 
mathematical framework, there are a large number of successful attempts in the literature to by-pass the limitations 
of perturbation theory. In particular, there are methods to perform resummations of perturbative series (even if 
they are divergent) , which ammounts in some cases to analytically continue weak-coupling series to a strong-coupling 
domain [12-19]. For instance, starting from a function of a coupling constant g defined formally by means of a scries 
(not necessarily convergent), 

oo 
n=0 

we can, under certain analyticity assumptions (the validity of the Watson-Nevanlinna-Sokal theorem, see for instance 
[20] and other references therein), define its Borel transform as the associated series, 

oo 

B(b) =]T ^6", (1.2) 

n=0 

which has an analytic continuation on a strip along the real 6-axis from zero to infinity. It can be easily verified that 
the inverse Borel transform 

1 f°° 

B(g) = - db e b /°B(b) , (1.3) 
9 Jo 

reproduces formally the original series (1.1). From a physical point of view, the important remark is that the series 
B(b) can be convergent and summed up even if the series (1.1) diverges. In this case, the inverse Borel tranform (1.3) 
defines a function of g, B(g), which we can think about as the "sum" of the divergent series (1.1). This function 
B{g) can be defined for values of g not necessarily small and, in this sense we can perform an analytic continuation 
from a weak to a strong-coupling regime. Techniques of this type are of a predictive character and have been largely 
employed in the last years in Quantum Field Theory literature. 

Nevertheless as a matter of principle, due to the non vanishing of the coupling constant, the idea of a bare particle 
associated to a bare matter field is actually an artifact of perturbation theory and is physically meaningless. A charged 
physical particle is always coupled to the gauge field, in other words, it is always "dressed" by a cloud of quanta of 
the gauge field (photons, in the case of Electromagnetic field). In fact as mentioned above, from a phenomenological 
point of view there are situations even in the scope of QED, where perturbation methods are of little usefulness, for 
instance, resonant effects associated to the coupling of atoms with strong radiofrcquency fields. As remarked in [3], 
the theoretical understanding of these effects using perturbative methods requires the calculation of very high-order 
terms in perturbation theory, what makes the standard Feynman diagrams technique practically unreliable. The trials 
of treating non-perturbativelly systems of this type, have lead to the idea of " dressed atom" , introduced originally in 
Refs. [21] and [22]. Since then this concept has been used to investigate several situations involving the interaction of 
atoms and electromagnetic fields, as for instance, atoms embedded in a strong radiofrequency field background and 
atoms in intense resonant laser beans. 

In order to give a precise mathematical definition and a clear physical meaning to the idea of dressed atom, we do 
not face the non-linear character of the problem involved in realistic situations, we consider instead, as mentioned 
above a linear problem, in which case a rigorous definition of "dressed atom" and more generally of a dressed particle 
interacting weakly or strongly with the field can be given. In other words, in this paper we adopt a general physicist's 
point of view. We do not intend to describe all the specific features of a real non-linear physical situation. Instead 
we analyse a simplified linear version of the atom-field or particle-environment system and we try to extract the most 
detailed information we can from this model. We will introduce dressed states, by means of a precise and rigorous 
definition to solve our problem. Our dressed states can be viewed as a rigorous version of the semiqualitative idea of 
dressed atom mentioned above, which can be constructed in reason of the linear character of our problem. Our aim 
in taking a linear model is to have a clearer understanding of what is one of the essential points, namely, the need for 
non-perturbative analytical treatments of coupled systems, which is the basic problem underlying the idea of a dressed 
quantum mechanical system. Of course, such an approach to a realistic non-linear system is an extremely hard task, 
and here we achieve what we think is a good agreement between physical reality and mathematical reliability. 



In recent publications [23,24] a method (dressed coordinates and dressed states) has been introduced that allows a 
non-perturbative approach to situations of the type described above, provided that the interaction between the parts 
of the system can be approximated by a linear coupling. More precisely, the method applies for all systems that can 
be described by an Hamiltonian of the form, 



H 



N 



fe=i 



N 



^2c k q k , 



(1.4) 



fc=i 



where the subscript refers to the particle and k = 1,2, ...N refer to the harmonic environment modes. The limit 
TV — > oo in Eq.(1.4) is understood. The equations of motion are, 



q (t) + UoQo(t) = c k1k(t) 



k=i 



and 



q k (t) + u k q k (t) = c k q (t) 



(1.5) 



(1.6) 



A Hamiltonian of the type above, describing a linear coupling of a particle with an environment, has been used 
in Refs. [8,9,25-28] to study the quantum Brownian motion of a particle. In the case of the coupled atom field 
system, this formalism recovers the experimental observation that excited states of atoms (weakly) coupled to the 
electromagnetic field in sufficiently small cavities are stable. It allows to give formulas for the probability of an atom 
to remain excited for an infinitely long time, provided it is placed in a cavity of appropriate size [24]. For an emission 
frequency in the visible red, the size of such cavity is in agreement with experimental observations [30,31]. The model 
that we consider for the particle-field system consists of an one dimensional oscillator interacting with a massless 
scalar field and described by the Lagrangean: 
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2ny/g j d 3 x(j>(x)6(x)q , 



(1.7) 



where g is a coupling constant with dimension of frequency 

Notice that we are supposing that the particle interacts with the field by a contact term only at the origin or, in 
other words the particle is centered at x = 0. This is equivalent to the dipole approximation for electromagnetic 
interactions. Next we confine the whole system into a sphere of diameter L. Writing 



4>{x, t) = ^2q k (t)(t> k (x) 



(1.8) 



where the 4> k (x) form an orthonormal basis obeying the equation, 

-V 2 fe (.f) = uj 2 k 4> k {x) 
replacing Eq. (1.8) in Eq. (1.7) and using Eq. (1.9) we obtain 



C = 



1 



•2 2 2 i W -2 2 2\ 



+ 2TTy/gq ^2^k(0)qk 



(1.9) 



(1.10) 



Notice that the field modes that interact with the particle arc those modes such that 0(0) ^ 0. These modes, 
given by the eigenvalues equation (1.9), possess spherical symmetry. Then, solving the eigenvalues equation (1.9) in 
spherical coordinates, with the boundary condition <f> k (L/2) = 0, we get for the spherically symmetric modes 
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sin(o;fer) 
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T~ i fe — 1, 2, 3, 



(1.11) 



where r — \x\. From the Lagrangean (1.10) we can get the Hamiltonian in a standard way. Using Eq. (1.11) we 
obtain the Hamiltonian given by Eq. (1.4) with the coupling coefficients c k — rju k , where r\ = \j2gKui and Aui is the 
interval between two neighbouring bath (field) frequencies. 



II. THE EIGENFREQUENCIES SPECTRUM AND THE DIAGONALIZING MATRIX 



The bilinear Hamiltonian (1.4) can be turned to principal axis by means of a point transformation, 
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«M = IXQr, ?V = E^' M=(0, {*:}), k = l,2,...,N 



(2.1) 



r=0 



performed by an orthonormal matrix T = (ip. The subscripts \i = and fj, = k refer respectively to the particle and 
to the harmonic modes of the bath and r refers to the normal modes. In terms of normal momenta and coordinates, 
the transformed Hamiltonian in principal axis reads, 



1 N 



(2.2) 



r=0 



where the fi r 's are the normal frequencies corresponding to the possible collective stable oscillation modes of the 
coupled system. Using the coordinate transformation = t^Q r in the equations of motion and explicitly making use 

of the normalization condition X^o(*p 2 = 1> we § ct > 
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(2.3) 



and the normal frequencies O r are given as solutions of the equation, 
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Remembering Ck — r\^k, Eq. (2.4) can be written as 



c 2 - Nrf -n 2 r = V 2 J2 



2 



(2.4) 



(2.5) 



In the limit TV — > oo Eq.(2.5) is meaningless if uj is finite. To overcome this difficulty we define ujq as containing 
the divergent term Nrf and a finite part, 



ujI = uj 2 + Nrj 2 



(2.6) 



where ui 2 is finite and is defined as the physical (squared) frequency, while uu 2 is not physical, it as a bare (squared) 
frequency and defined in such a way that the divergent terms in the left hand side of Eq. (2.5) cancel. Replacing 
Eq.(2.6) in Eq. (2.5) we obtain 



N 



-n 2 r = v 2 J2 



n 2 



(2.7) 



We see that the above procedure is exactly the analogous of naive mass renormalization in quantum field theory: the 
addition of the counterterm —(r) 2 N/2)q 2 allows to compensate the infinity of to 2 , in such a way as to leave a finite, 
physically mcaninful renormalized frequency uj. This simple renormalization scheme has been originally introduced 
in Rcf . [32] . For a nice discussion on the subjec the interested reader can see Ref. [33] . Using the formula, 



y 1 



Eq. (2.7) can be written in closed form, 



,LO n 2 , Ltd 2 , 



(2.8) 



(2.9) 



The solutions of Eq. (2.9) with respect to Q give the spectrum of eigenfrequencies Sl r corresponding to the collective 
normal modes. The transformation matrix elements T = (t r ) are obtained in terms of the physically meaningful 
quantities O r , u>, after some rather long but straightforward manipulations. They read, 

t r - ^ r _ rpjk r , , 

sJ(Q? r - (D 2 ) 2 + ^(30 2 - CJ 2 ) + TrVft 2 W * - "r 

In free-space, that is, in the limit L^oowe get for t^, 



and for t£, 



= lim - v y 2.11) 

AO^O ^(fi 2 - w 2 ) 2 + ttV^ 2 



t r = 2gu k Au tl r 

k {Uj k + fl r )(0J k - Sir) ^/(O 2 - UJ 2 f + 7T 2 . 9 2 2 ' 



where we use the fact that in the limit L — > 00, AO — > Aoj = 2ir/L. 



III. DRESSED STATES AND THE EMISSION PROCESS IN FREE SPACE 



We define below some coordinates q' , q[ associated respectively to the dressed particle (atom) and to the dressed 
bath (field). These coordinates will reveal themselves to be appropriate to give an appealing non-perturbative de- 
scription of the particle-field system. The normalized eigenstates of our system (eigenstates of the Hamiltonian in 
principal axis) can be written in terms of the normal coordinates, 



(<9|n ,ni, ...;t) 




r (QK 



(3.1) 



s=0 



where H ns stands for the n s -th Hermite polynomial and T {Q) is the normalized ground state eigenfunction. The 
eigenvalues E noni ... are given by 



En 
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E 

s=0 



n. 



(3.2) 



Next we intend to divide the system into the dressed particle and the dressed environment (bath, field) by means 
of some conveniently chosen dressed coordinates, q' and q'j associated respectively to the dressed particle and to 
the dressed oscillators composing the environment. These coordinates will allow a natural division of the system as 
composed of the dressed (physically observed) particle and the dressed environment. The dressed particle will contain 
automatically all the effects of the environment on it. Clearly, these dressed coordinates should not be introduced 
arbitrarilly. Since our problem is linear, we will require a linear transformation between the normal and dressed 
coordinates (different from the transformation (2.1) linking the normal to the bare coordinates). Also, we demand 
the physical condition of vacuum stability. We assume that at some given time (t = 0) the system is described by 
dressed states, whose wavefunctions are defined by, 
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vw-.^o) = n 
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W) 



(3.3) 



where q'^ = (q'o,q'i), to^ — (u),u>i) and To(q') has the same functional dependence on the variables q' as the ground 
state eigenfunction r (Q) has on the variables Q, i.e., can be obtained from T (Q) by replacing Q r and f2 r by the 
dressed variables q'^ and u)^. The dressed states given by Eq. (3.3) describe the dressed particle in its fco-th excited 
level in the precensse of h dressed field quanta (photons in the case of electromagnetic interactions) of frequencies 
u>i. Note that the above wavefunctions will evolve in time in a more complicated form than the unitary evolution of 
the eigenstates (3.1), since these wavefunctions are not eigenstates of the diagonal Hamiltonian (1.4). 



In order to satisfy the physical condition of vacuum stability (invariance under a tranformation from normal to 
dressed coordinates) we remember that the the ground state eigenfunction of the system has the form, 

r (Q) oc e'^ Sr=o n <-Qr ; (3.4) 
and we require that the ground state in terms of the dressed coordinates should have the form 

roteOoce^^o"'^ 2 . (3.5) 

From Eqs. (3.4) and (3.5) it can be seen that the vacuum invariance requirement is satisfied if we define the dressed 
coordinates by, 

AT 

v^g^IXv^VQr. (3.6) 

These dressed coordinates are new collective coordinates, different from the bare coordinates describing the bare 
particle and the free field modes, and also from the normal (collective) coordinates {Q r }- Indeed, using Eq. (2.1) in 
Eq. (3.6) the dressed coordinates can be written in terms of the bare coordinates as 



N N 

r = y2ai llv q v , a^v = —= Y^CV^V (3.7) 



r=0 



Our dressed states, given by Eq.(3.3), are collective but non stable states (only the the lowest energy state, the 
invariant vacuum, is stable), linear combinations of the (stable) eigensatates (3.1) defined in terms of the normal 
coordinates. The coefficients of these combinations are given in Eq. (3.9) below and explicit formulas for these 
coefficients for an interesting physical situation are given in Eq.(3.12). This gives a complete and rigorous definition 
of our dressed states. Moreover, our dressed states have the interesting property of distributing the energy initially in a 
particular dressed state, among itself and all other dressed states with precise and well defined probability amplitudes 
[23]. We choose these dressed states as physically meaningful and we test successfully this hypothesis by studying 
the radiation process by an atom in a cavity. In both cases, of a very large or a very small cavity, our results are in 
agreement with experimental observations. 

Using Eq. (3.6) the functions given by Eq. (3.3) can be expressed in terms of the normal coordinates Q r , but since 
the functions <p no ,n 1 ....(Q) form a complete set of orthonormal functions, the functions V'fco,fci,...(<7 / ) can be written as 
linear combinations of the eigenfunctions of the coupled system (we take t — for the moment), 

oo 

Vw.(<z',o) = E T "o°" 1 1 .::>«o™ 1 ...(Q,o), (3.8) 

noni ...—0 

where the coefficients are given by, 

TZZ 1 :;; = J dQ^ KoK1 ...( q ')<Pn oni ...(Q) , (3-9) 

the integral extending over the whole Q-space. We consider the particular configuration tp in which only the dressed 
particle (atom) q' is in its fc-th excited state, all other being in the ground state, 

Vfeoo...(g',0) = (2- k k\)-iH k (VQq' )T (q') . (3.10) 

The above initial state represents a particle in the fc-th excited level and no field quanta. The coefficients given by 
Eq. (3.9) can be calculated in this case using the theorem [34], 
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Hk ( J__ — _ r . ) = J- ^ t | ••• g mo (7^Q )g mi (v / »iQi)-. (3.11) 

y Z^rV V' m a +m 1 +... = k 



Replacing Eq. (3.10) in Eq. (3.9) and using the theorem given by Eq. (3.11), we get 



Now we would like to compute the probability amplitude that at time t the atom still remain in the A:-th excited level, 
that we denote by Ak(t), and given by 

A k (t) = (Vfcoo..(g',0), Vfcoo... (</,*)) , (3.13) 

where Vfeoo...(9'j t) can be obtained from Eq. (3.8), 

oo 

V>fcoo... (q',t)= J2 T^--e-* E ^-^ non UQ,0)- (3-14) 

noni ...—0 

Replacing Eq. (3.14) in Eq. (3.13) we get 

M*)= E ( T E-) 2e " iVl - ( 3 - 15 ) 

no,ni,...=0°° 

Taking k = 1 in Eq. (3.12) and replacing in Eq. (3.15), we obtain for the probability amplitude that at time t the 
atom still remain in the first excited level, 

N 

f(^E(« 2 ^ rt ' ( 3 - 16 ) 

r=0 

where we do not include a global phase factor that will not contribuite to the probability |/ 00 (i)| 2 and we have changed 
the notation f 00 (t) = Ai(t). In the case of a very large cavity (L — > oo) the probability amplitude that the dressed 
particle be still excited at time t can be obtained using Eq. (2.11) to transform the discrete sum given by Eq. (3.16) 
into an integral, 



/•oo q2 —iQt 

r(t)=^J o dn — U Z ; (3.1.7) 



'(n 2 -oj 2 ) 2 + n 2 g 2 n 2 ■ 



For large t (t >> i), but for in principle arbitrary coupling g, we obtain for the probability of finding the atom still 
excited at time t, the result [35], 



16g 2 

+ -4' (3 - 18) 



where k — \J uj 2 — ^f-. In the above expression the approximation t >> i plays a role only in the two last terms, due 

to the difficulties to evaluate exactly the integral in Eq. (3.17) along the imaginary axis using Cauchy's theorem. The 
first term comes from the residue at O = k + i^- and would be the same if we have done an exact calculation. If we 
consider in eq. (3.18) g « u>, which corresponds in electromagnetic theory to the fact that the fine structure constant 
a is small compared to unity (for explicit calculations we take below g — aui), we obtain the well known perturbative 
exponential decay law. For arbitrary time and for arbitrary ratio g/ui we can integrate Eq. (3.17) numerically and 
we obtain an exponential like decay law. 

Remark that in the limit of a very large cavity, we can express the dressed coordinates q'^ in terms of the bare ones, 

in the limit L — > oo, 

q'o = A 00 (u),g)qo, (3.19) 
q'i = Qi, (3.20) 



where Aoo(oj,g), is given by. 



1 f 00 2gn 2 VUdQ 
Aoo^g),-^^ (02 __ 2)2+7rW - (3.21) 



It is interesting to compare Eqs.(3.7) with Eqs.(3.19), ( 3.20). In the case of Eqs.(3.7) for finite R, the coordinates 
q' and {q-} are all dressed, in the sense that they are all collective, both the field modes and the particle can not be 
separated in this language. In the limit R — ► oo, Eqs.(3.19) and (3.20) tells us that the coordinate q' describes the 
particle modified by the presence of the field in a indissoluble way, the particle is always dressed by the field. On 
the other side, the dressed harmonic modes of the field, described by the coordinates q\ are identical to the bare field 
modes, in other words, the field keeps in the limit R — > oo its proper identity, while the particle is always accompanied 
by a cloud of field quanta. Therefore we identify the coordinate q' as the coordinate describing the particle dressed 
by its proper field, being the whole system divided into dressed particle and field, without appeal to the concept of 
interaction between them, the interaction being absorbed in the dressing cloud of the particle. In the next section we 
will consider a partcle (an atom) confined in a small cavity. By "small" it is understood a cavity having dimensions 
much lower than macroscopic dimensions, but still much larger than atomic dimensions. In this case the concepts 
of dressed atom and of dressed field can be taken as the physically meaningful ones to describe the atom radiation 
process [29]. 



IV. THE CONFINED SYSTEM 



Let us now consider the particle (an atom) placed in the center of the cavity in the case of a small diameter L, 
i.e. that satisfies the condition of being much smaller than the coherence lenght, L << 2c/ g [29]. To obtain the 
eigenfrequencies spectrum, we remark that from an analysis of Eq. (2.9) it can be seen that in the case of a small 
values of L, its solutions are near the frequency values corresponding to the asymptots of the curve cot(^-), which 
correspond to the field modes Wj = 2ni/L. The smallest solution departs more from the first asymptot than the other 
larger solutions depart from their respective nearest asymptot. As we take larger and larger solutions, they are nearer 
and nearer to the values corresponding to the asymptots. For instance, for a value of L of the order of 2 x 10~ 2 m 
and u> ~ 10 10 /s, only the lowest eigenfrcquency Oo is significantly different from the field frequency corresponding 
to the first asymptot, all the other eigenfrequencies fife, k = 1,2, ... being very close to the field modes 2wk/L. For 
higher values of u> (and lower values of L) the differences between the eigenfrequencies and the field modes frequencies 
are still smaller. Thus to solve Eq. (2.9) for the larger eigenfrequencies we expand the function cot(4p) around the 
values corresponding to the asymptots. Writing, 

2ir 

n k = —(k + e k ), fc = l,2,.. (4.1) 
with < efe < 1, and substituting in Eq. (2.9) we get, 

cot(7r6 fe )-4(fc + £fe ) + -^(l-^). (4.2) 
gL (k + e k ) 2irg 

But since for a small value of L every e k is much smaller than 1, Eq. (4.2) can be linearized in e kl giving, 

_ AvgLk 

€k 2(4^fc 2 -^ 2 £ 2 )- { ' 

Eqs. (4.1) and (4.3) give approximate solutions for the eigenfrequencies fl k , k = 1, 2.... To solve Eq. (2.9) with respect 
to the lowest eigenfrequency flo, let us assume that it satisfies the condition floL/2 << 1 (we will see below that this 
condition is compatible with the condition of a small L as defined above). Inserting the condition Q.qL/2 << 1 in Eq. 
(2.9) and keeping up to quadratic terms in O we obtain the solution for the lowest eigenfrcquency flo, 

O = ; = (4.4) 

2c 

Consistency between Eq. (4.4) and the condition fl L/2 « 1 gives a condition on L, 

2c7T 



9 2 



1 + Wi + 



ng 



(4.5) 



a) Weak coupling 



Let us particularize our model to situations where interactions of electromagnetic type are involved. In this case, 
let us define the coupling constant g to be such that g = Qa, where a is the fine structure constant, a = 1/137. 
Then the factor multiplying 2/g Eq. (4.5) is <~ 0.07 and the condition L <C 2/g is replaced by a more restrictive one, 
L <C 0.07(2/g). For a typical infrared frequency, for instance w~2,0x 10 n /s, our calculations are valid for a value 
of L, L <C 10~ 3 to [29]. From Eqs.(2.10) and using the above expressions for the eigenfrequencies for small L, we 
obtain the matrix elements, 

«8) a «i-^, (*S) a «S- (4 - 6) 

To obtain the above equations we have neglected the corrective term e k , from the expressions for the eigenfrequencies 
ilk- Nevertheless, corrections in e k should be included in the expressions for the matrix elements t%, in order to avoid 
spurious singularities due to our approximation. Let us consider the situation where the dressed atom is initially in 
its first excited level. Then from Eq. (3.16) we obtain the probability that it will still be excited after a ellapsed time 

t, 

OO OO 

l/°°(*)| 2 - (O 4 + 2£(*o) 2 (*o) 2 cos(n k - n )t + J2 (*o) 2 (*o) 2 cos(tt fc - fi,)t , (4.7) 

fe=l k,l=l 

or, using Eqs.(4.6) in Eq.(4.7), 

|/ 00 (i)| 2 « 1 - irS + 4(£ - 5 2 ) £ ^ cos(n k - Qo)t + n 2 5 2 + -±<5 2 ^ cos(n k - Sl t )t, (4.8) 

fc=l k,l=l 

where we have introduced the dimensionless parameter 6 = Lg/2 -C 1, corresponding to a small value of L and we 
remember that the eigenfrequencies arc given by Eqs. (4.1) and (4.3). As time goes on, the probability that the 
mechanical oscillator be excited attains periodically a minimum value which has a lower bound given by, 

Min(|f 00 (t)| 2 ) = l-f* + ^ 2 . (4.9) 

For a frequency uj of the order u> <~ 4.00 x 10 14 /s (in the red visible), which corresponds to 5 ~ 0.005 and L <~ 
1.0 x 10~ 6 m [29], we see from Eq.(4.9) that the probability that the atom be at any time excited will never fall 
below a value <~ 0.97, or a decay probability that is never higher that a value ~ 0.03. It is interesting to compare 
this result with experimental observations in [30,31], where stability is found for atoms emiting in the visible range 
placed between two parallel mirrors a distance L = 1.1 x 10~ 6 m apart from one another. For lower frequencies the 
value of the spacing L ensuring quasi-stability of the same order as above for the excited atom, may be considerably 
larger. For instance, for u in a typical microwave value, ui ~ 2, 00 x 10 10 /s and taking also S ~ 0.005, the probability 
that the atom remain in the first excited level at any time would be larger than a value of the order of 97%, for 
a value of L, L <~ 2.0 x 10~ 2 m. The probability that the mechanical oscillator remain excited as time goes on, 
oscillates with time between a maximum and a minimum values and never departs significantly from the situation 
of stability in the excited state. Indeed for an emission frequency Co ~ 4.00 x 10 14 /s (in the red visible) considered 
above and L ~ 1.0 x 10 _6 m, the period of oscillation between the minimum and maximum values of the probability 
that the dressed mechanical oscillator be excited, is of the order of x 10~ 14 s, while for Q ~ 2, 00 x 10 10 /s, and 
L~2.8x 10~ 2 m, the period is of the order of ^ x 10~ 10 s. 
b) Strong coupling 

In this case it can be seen from Eqs (4.3), (4.4) and (4.5) that £1 ~ w and we obtain from Eq. (2.10), 

< ! S> 2 *ttW? (4 - 10) 

Using Eqs. (4. 10) in Eq.(4.7), we obtain for the probability that the excited particle be still at the first energy level at 
time t, the expression, 

(2 \ 2 9 °° ox a °° i 

2TVs) + —s E w cos(0fe - ^ + ^ 2 £ m cos( ^ - Ql)t • (4 - n) 
7 k=i k.i=i 



We see from (4.11) that as the system evolves in time, the probability that the particle still be excited, attains 
periodically a minimum value which has a lower bound given by, 



The behaviour of the system in the strong coupling regime is completely different from weak coupling behaviour. 
The condition of positivity of (4.12) imposes for fixed values of g and Q an upper bound for the quantity 5, S max , 
which corresponds to an upper bound to the diameter L of the cavity, L max (remember S = Lg/2). Values of S larger 
than 5 max , or equivalently, values of L larger than L max are unphysical (correspond to negative probabilities) and 
should not be considered. These upper bounds are obtained from the solution with respect to 5, of the inequality, 
Min(<5) > 0. We have Min(5) > or Min((5) = 0, for respectively S < 8 max or S — 5 max . The solution of the equation 

Min((5) = gives S max = — ^ + V~ 2 +3\/5 w Q 37 p or § _ S max , the minimal probability that the excited particle 
remain stable vanishes. We see comparing with the results of the previous subsection, that the behaviour of the 
system for strong coupling is rather different from its behaviour in the weak coupling regime. For appropriate cavity 
dimensions, which are of the same order of those ensuring stability in the weak coupling regime, we ensure for strong 
coupling the complete decay of the particle to the ground state in a small ellapsed time. 



V. AN ANHARMONIC OSCILLATOR IN A LARGE CAVITY 



In this section we intend to generalize, making some appropriate approximations,to an anharmonic oscillator, the 
linear coupling to an environnement as it has been considered in the preceding sections. For more details see Ref. 
[39]. In this case, the whole system is described by the Hamiltonian, 

N 

H[X r ] = H(p , q , {p k , q k }) + E KT^q^q^ , (5.1) 

r=0 

where H(j>q, go, {pk, qk}) is the bilinear Hamiltonian given by Eq. (1.4) and T p l) pa are some coefficients that will be 
defined below. In Eq. (5.1) summation over repeated greek labels is understood. We do not intend to go to higher 
orders in the perturbative series for the energy eigenstatates, we will remain at a first order correction in A r , and we 
will try to see what are the effects of the anharmonicity term, given by the last term of Eq. (5.1), on our previous 
results for the linear coupling with an environnement (field). We notice that the anharmonicity term in Eq. (5.1) 
involves, independent quartic terms of the type /1 = 0, {i} (self-coupling of the bare oscillator and of the field 
modes), quartic terms coupling the oscillator to the field modes and also the terms coupling the field modes among 
themselves. These terms are of the type qoqf, 9o*?f > QoQi' an d of the type qfqj 1 , with n + m = 4, for the coupling 
between the field modes. We intend to start from the exact solutions we have found in the linear case, and investigate 
at first order, how the quartic interaction characteristic of the anharmonicity changes the decay probabilities obtained 
in the linear case. 

The introduction of the quartic interaction term in Eq.(5.1) changes the Hamiltonian in principal axis from Eq. 
(2.2) into 

N N 

h[K] = ^ E( p ' + n rQr) + E Kr^%n r ;QQ ri Q r2 Q r3 Q r4 , (5.2) 

where summation over the repeated indices V\,ri,r%,r\ is understood. In order to have a specific quartic interaction, 

(r) 

we make a choice for the coefficients T^ vpcr in the above equation, 

q-(r) _ f r f r f r f r /c o\ 

which replaced in Eq. (5.2) and using the orthonormality of the matrix t 1 ^ gives the Hamiltonian in principal axis, 

N N 

H [A r ] - - + ^0?) + E X rQt ■ ( 5 - 4 ) 

r=0 r=0 

Performing a perturbative calculation in A r , we can obtain the first order correction to the energy of the state 
corresponding to n — n\ = .... = 1 in Eq. (3.2), in such a way that the A r -corrected energy of such state can be 
written in the form, 



JV 

E[X r ] = + Ar4 X) ) , 

r=0 



(5.5) 



where we do not include the energy of the ground state for the reason explained after Eq. (3.16). The first order 
correction in Eq. (5.5) is given by 

4" - ^ • (5-6) 

For sufficiently small A r we will have quasi-harmonic normal collective modes having frequencies tt r + A r e[ 1 ' ) . Accord- 
ingly we can describe approximately the system in terms of modified harmonic cigenstates, which can be written as a 
generalization of the exact eigenstates (3.1), replacing the eigenfrequencies Q r by the A-correctcd values il r + \ r e"r\ 



N 



XQ,t;{X s }) = Yl 



s=0 



lon s n 
! —H ns {f-fQs) 



n.,! 



r (g)e~ l S 3 " s ( Os+Ase ' 1) ) t . (5.7) 



From the modified harmonic eigenstates (5.7), we can follow analogous steps as in the harmonic case to study the 
{A r }-corrected evolution of a dressed particle, generalizing Eq. (3.16), obtaining 



N 



/°°(t; {A r }) - ^(g2 e -(SWA r eW)t _ (5 8) 

For a very large cavity, L — > oo, using Eq. (2.11), we obtain we obtain the continuum limit, 

/°°(*;{An}) = 2,/ o dn — __— - . (5.9) 

From dimensional arguments, we can choose Aq = AS1 3 , where A is a dimensionless small fixed constant. With this 
choice, after expanding in powers of A the exponential in Eq. (5.9), we obtain to first order in A the amplitude, 

/°° ( i ; A) = />) + ^|/ 00 (t), (5.10) 

where f 00 (t) is the probability amplitude, for the harmonic problem, that the particle be after a time t, still in the 
first excited level. From Eq. (5.10) we can obtain at order A, the probability that the particle remains in the first 
excited state at time t, 

|/°° ( i;A)| 2 H/>)| 2 + ^|l/°W- (5.11) 

We know that |/ 00 (i)| 2 is a decreasing function of t, what means that the derivative of this function with respect 
to t is negative. Therefore, since A > 0, we conclude from Eq. (5.11) that |/ 00 (i;A)| 2 is smaller than the harmonic 
probability |/ 00 (i)| 2 - In Fig. 1 we plot on the same scale the A-corrected probability (5.11) and the harmonic 
probability given by Eq. (3.18), for u> = 4.0 x 10 14 /s and g = aoj, where a is the fine structure constant, a = 1/137 
and the time is rescaled as t x 10 _13 s. The solid line is the harmonic probability (3.18) and the dashed line is the 
A-corrected probability (5.11), for A = 1/50. We see clearly the enhancement of the particle decay induced by the 
quartic interaction. 
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FIG. 1. Plot on the vertical axis, commonly named P[t], of the A-corrected probability (5.11) (dashed line) and the harmonic 
probability Eq.(3.18) (solid line), for Q — 4.0 x 10 14 /s and g = auj. a is the fine structure constant, a = 1/137 and time in 
units of 10" 13 s 



VI. CONCLUSIONS 

In this paper we have analysed a linearized version of an particle-environment system and we have tried to give the 
more exact and rigorous treatment we could to the problem. We have adopted a general physicist's point of view, in 
the sense that we have renounced to approach very closely to the real behaviour of a complicated non-linear system, 
as a quark-gluons coupled system, to study instead a linear model. As a counterpart, an exact solution has been 
possible. 

We have presented a summary of mathematical results previously obtained, to describe an ohmic quantum system 
consisting of a particle (in the larger sense of a " material body" , an atom or a Brownian particle for instance) coupled 
to an environment modeled by non-interacting oscillators. We have used a formalism (dressed coordinates and dressed 
states) that allows a non-perturbative approach to the time evolution of the system, in rather different situations as 
confinement or in free space. For finite R, our system could be seen as a simplified linear model to confined quarks 
and gluons inside a hadron. In this case all dressed coordinates are effectively dressed, in the sense that they are 
all collective, both the field modes and the particle can not be separated in this language [see Eq.(3.7)]. Of course 
the normal coordinates are also collective, but they correspond to stable eigenstates, no change in time exists for 
them. If we ascribe physical meaning to our dressed coordinates and states in the context of our model, matter and 
gauge quanta inside hadrons, can not be individualized as " quarks" and " gluons" , all we have is a kind of quark-gluon 
"magma". Since quarks and gluons are permanently confined, we could roughly think that, in the context of the 
model studied here, quarks and gluons should not really exist, they would be an artifact of perturbation theory. 

In the limit R — > oo, we get the result that the dressed coordinate associated to the particle describes the particle 
modified by the presence of the field in a indissoluble way, the particle is always dressed by the field. On the other 
side, the dressed harmonic modes of the field, are in the limit R — > oo identical to the bare field modes, in other 
words, the field keeps in the free space limit its own identity, while the particle is always accompanied by a cloud of 
field quanta [see Eqs.(3.19) and (3.20)]. 

The study on the behaviour of particles (for instance atoms in the harmonic approximation) confined in small 
cavities, shows that it is complctelly different from the behaviour in free space. We have implicitly assumed in this 
study that a small cavity is still much larger than atomic dimensions (which is indeed the case for the experimental 
situation compared to our results, corresponding to a cavity diameter of ~ in such a way that the dressed particle 
could be a good approximation to the atom inside the cavity. In the first case the time evolution is very sensitive to 
the presence of boundaries, a fact that has been pointed out since a long time ago in the literature ( [36] , [37] , [38] ) . 
Our dressed states approach gives an unified description for the dressing of a charged particle by the field modes and 
the time evolution in a cavity of arbitrary size, which includes microcavities and very large cavities (free space). If 
we assume that our dressed particle is a good representation for an atom under certain circumstances, we recover 
here with our formalism the experimental observation that excited states of atoms in sufficiently small cavities are 
stable for weak coupling. In the weak coupling regime, we are able to give formulas for the probability of an atom to 
remain excited for an infinitely long time, provided it is placed in a cavity of appropriate size. For an interaction of 
electromagnetic type, for an emission frequency in the visible red, the size of such cavity is in good agreement with 
experimental observations. Also, our approach gives an exact result for emission in free space, generalizing the well 



known exponential decay law. The behaviour of the system for strong coupling is rather different from its behaviour in 
the weak coupling regime. For appropriate cavity dimensions, which are of the same order of those ensuring stability 
in the weak coupling regime, we ensure for strong coupling the complete decay of the particle to the ground state in a 
small ellapsed time. One possible conclusion is that by changing conveniently the physical and geometric parameters 
(the emission frequency, the strength of the coupling and the size of the confining cavity) our formalism theoretically 
allows a control on the rate of emission and of the energy storage capacity (perhaps information?) in the cavity. 
Depending on the strength of the coupling, the emission probability ranges from a complete stability to a very rapid 
decay. 

The formalism that we presented is rather general an applies to all systems that can be modeled by a Hamiltonian 
of the type given by Eq. (1.4). Applications to the termalization and decoherence problems will be reported elsewhere 
[40,41]. 
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